The formation of carotid artery plaque is closely related to complex hemodynamic factors. The accurate simulation of complex flow conditions is of great significance for the clinical diagnosis of carotid artery plaque. In order to simulate the pulsating flow accurately, Large Eddy Simulation (LES) model was combined with Lattice Boltzmann Method (LBM) to construct a LBM-LES carotid artery simulation algorithm, and a real geometric model of carotid artery stenosis was established through medical image reconstruction software, thus the high-resolution numerical simulation of carotid artery stenosis flows was conducted. By calculating blood flow velocity and Wall Shear Stress (WSS), some meaningful flow results were obtained, proving the effectiveness of LBM-LES in the study of blood flow in the carotid artery narrow posterior. Based on the OpenMP programming environment, the parallel computation of the grid of ten million magnitude was carried out on the fully interconnected fat node of high-performance cluster machine. The results show that the LBM-LES carotid artery simulation algorithm has good parallel performance.